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9^ ABSTRACT 

■ The luminosity distance vs. redshift law is now measured using supernovae and 7- 

ray bursts, and the angular size distance is measured at the surface of last scattering 
j> ■ by the CMB and at z = 0.35 by baryon acoustic oscillations. In this paper this data 

is fit to models for the equation of state with w = —1, w = constant, and w(z) = 
wo + w a (l — a). The last model is poorly constrained by the distance data, leading to 
unphysical solutions where the dark energy dominates at early times unless the large 
scale structure and acoustic scale constraints are modified to allow for early time dark 
energy effects. A flat ACDM model is consistent with all the data. 
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1. INTRODUCTION 

^ ■ Wood-Vasey et al. (2007) recently published supernova data from the ESSENCE project, while 

Riess et al. (2007) have published a large sample of supernovae from the SNLS project (Astier et al. 
2006), the HST, and Hi-z Supernova Team. Schaefer (2006) has published a sample of 7-ray burst 
distances. While GRBs give much less accurate distances than supernovae, they extend to much 
higher redshifts and the GRB data helps to distinguish between non-flat geometries and equations 
of state with i»/-l which both affect the distance-redshift law at 0(z s ) for low redshift. 

The analysis of Wood-Vasey et al. (2007) plotted contours of w and w a , based only on a subset 
of the supernova data used here and a prior on fi^f. These contours extended into the region where 
the dark energy dominated the density at the surface of last scattering or during nucleosynthesis. 
Li et al. (2006) added GRBs and large scale structure data to the SNe data but still show contours 
extending into the early dark energy domination region. Both Barger et al. (2006) and Riess et al. 
(2007) used w{z) laws in which the deviation of w from —1 was terminated for z > 1.8, beyond the 
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Table 1: Mean distance modulus relative to a Milne model for the 358 supernovae in Riess et al. 
(2007) Gold+Silver and Wood-Vasey et al. (2007). 
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redshift of the most distant supernova in the sample, but this is an arbitrary limit which would have 
to modified to allow for GRBs. Alam, Sahni & Starobinsky (2006) analyzed the SNe data using a 
quadratic polynomial in (1 + z) as the form for pDE- This form for py)e will never be dominant 
at early times since the matter density varies like (1 + z) 3 . Davis et al. (2007) have analyzed a 
subset of the combined supernovae, and have used an approximation to the CMB acoustic peak 
constraint that fails when dark energy dominates at high z. The part of parameter space where 
dark energy dominates at high z obviously should be excluded, and I show in this paper that 
appropriate modifications to the standard formulae for the acoustic scale, the F parameter, and 
Big Bang nucleosynthesis (BBNS) to allow for the possible importance dark energy at z ~ 10 9 
or z ~ 10 3 will lead to this exclusion automatically. The BBNS, acoustic scale and T limits are 
given in general forms involving pde(z), which can be used for any form of the equation of state. 
But the w = w + w a {l — a) form adopted by the Dark Energy Task Force report (Albrecht et al. 
2006) is used for all the specific plots in this paper. This form can be used for all redshifts without 
introducing a redshift cutoff which would be a new and otherwise unnecessary parameter. 
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Table 2: Mean distance modulus relative to a Milne model for 272 supernovae in Riess et al. (2007) 
Gold and Wood-Vasey et al. (2007). 
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2. OBSERVATIONS 
2.1. Supernovae 

The distance modulus vs. redshift data from Riess et al. (2007) were taken from the Web site 
provided by Riess. The distance moduli and redshifts for the ESSENCE supernovae were extracted 
from Table 9 in the Latex file for |astro-ph/0701041| (Wood-Vasey et al. 2007). Typically different 



groups analyze supernovae with different assumptions about the Hubble constant or equivalently 
the absolute magnitude M. of a canonical SN la with a nominal decay rate. In order to combine the 
new supernovae from ESSENCE with the Riess et al. sample, it was necessary to check the relative 
normalization of the two data sets using the 93 objects they have in common. Figure [T] shows the 
comparison. The scatter in the differential distance moduli is 0.2 mag la, which seems unusually 
high, and the median difference in [i is 0.022 mag which is consistent with the standard deviation 
of the mean given the scatter. Objects which are not in the Riess et al. sample but which had 
successful fits with x 2 per degree of freedom < 7 were added to the Riess et al. sample, with the 
0.022 mag added to [A, and an intrinsic scatter of 0.10 mag added in quadrature to a^. This gives 
a total sample of 358 SNe. I have binned the SNe into bins containing < 2 + Ntot/W objects and 
widths less than 0.1 in redshift. An empty Universe model (the Milne model) was first subtracted 
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Table 3: Mean distance modulus relative to a Milne model for 69 GRBs in Schaefer (2007). 
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from the /u values. The binned values are listed in Table[TJ A Hubble constant of 63.8 km/sec/Mpc 
was used when computing the Milne model, but this value has no effect on the parameter limits 
computed in this paper. Its only effect is to add a constant to the {Afi) values in the Tables. The 
mean differential distance modulus in each bin is found by minimizing a modified 

X 2 = Yl f((< A ^ >~lf i i~ HMilne])/<?i) (1) 
i 

where f(x) = x 2 for \x\ < 2, or 4|x| — 4 otherwise. The modification deweights extreme outliers. 
Riess et al. (2007) recommend dropping SNe with redshifts less than 0.023 to avoid a possible 
"Hubble bubble" seen by Jha et al. (2007), but I have instead used a large velocity error of a v = 1500 
km/sec which gives an extra of (5/ In W)(a v /cz) which is added in quadrature with the tabulated 

Table [2] was constructed the same way but omitting the "Silver" objects in Riess et al. (2007). 

The binned values from Tables Q] and [2] are shown in Figure [2] along with the flat ACDM model 
that best fits the Hubble diagram data alone. This model has £Im = 0.369. There is an excursion 
around z = 0.5 that can be seen clearly in the binned supernova data. A simple 3 parameter fit to 
this bump gives a Ax 2 of 15 in the total sample, but only 6 if the "Silver" SNe are excluded. 

2.2. GRBs 

Schaefer (2007) has given a sample of 69 GRBs with redshifts and distance moduli. These 
values have been binned as well, but with bin widths < 0.1(1 + z m i n + z max ). The binned values 
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Fig. 1. — Difference in distance moduli for objects common to Wood-Vasey et al. (2007) (he) an d 
Riess et al. (2007) (hr). The RMS scatter is 0.2 mag. 



are listed in Table [3l In constructing the Table, a Milne model with a Hubble constant of 72 
km/sec/Mpc was subtracted from the individual distance moduli before minimizing the modified 
X 2 . Figure [3] shows the binned data from both the GRBs and the supernovae. 



2.3. Hubble Constant 

The Hubble constant largely cancels out in analyses of supernovae and 7-ray burst distances. 
But the Hubble constant does enter into converting the T parameter (£lMh) into a prior on Ojvf. 
The Hubble constant data used in this paper comes from Freedman et al. (2001), the DIRECT 
project double-lined eclipsing binary in M33 (Bonanos et al. 2006), Cepheids in the nuclear maser 
ring galaxy M106 (Macri et al. 2006), and the Sunyaev-Zeldovich effect (Bonamente et al. 2006). 
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Fig. 2. — Binned supernova data vs. redshift compared to a flat ACDM model with Oj^f = 0.369. 
The filled circles are binned points from the full dataset, while the open circles have omitted the 
"Silver" subset. 



These papers gave values of 72 ± 8, 61 ± 4, 74 ± 7 and 77 ± 10 km/sec/Mpc. Assuming that the 
uncertainties in these determinations are uncorrelated and equal to 10 km/sec/Mpc after allowing 
for systematics, the average value for H Q is 71 ± 5 km/sec/Mpc. This average is consistent with 
the 74 ± A{stat) ± 5(sys) km/sec/Mpc from Riess et al. (2005). 



2.4. Matter Density 

The matter density SIm^ 2 is fairly well determined by fitting the CMB power spectrum. In 
this paper the non-flat ACDM chain (ocdm_wmap_l.txt) at the LAMBDA data center has been 
used to determine average values for parameters. This chain gives com = ^l&ih 2 = 0.1289 ± 0.0079. 
The baryonic density is also well determined, with u>b = ^_b/i 2 = 0.02178 ± 0.00072. 
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Fig. 3. — Binned supernova and GRB data vs. redshift compared to a flat ACDM model with 
= 0.369, and non-flat ACDM model with Qm = 0.416 and Qtat = 1-115, and an exponentially 
evolving model in a £Im = 1 model (Wright 2002). 



2.5. Large Scale Structure 

Large scale structure data comes from the "big bend" in the power spectrum P(k). This has 
been measured by two large galaxy surveys: the SDSS and the 2dF. The SDSS gives a value for 
r = Q M h = 0.213 ± 0.0233 (Tegmark et al. 2004), while the 2dF gives (Cole et al. 2005) 

n M h = T true = 0.168 ± 0.016 + 0.3(1 - n s ) + 1.2fi„/0 M (2) 

I assume that the Tegmark et al. value uses n s = 1 and £l u = 0, and that it has the same sensitivity 
to these parameters as the 2dF. The combination of these two values gives a x 2 = 2.6 for 1 degree 
of freedom. While this is higher than the expected value of 1, it is certainly not high enough to 
trigger grave concerns. The neutrino density is uncertain but the minimal hierarchical mass pattern 
gives n v /n M » 0.004 while the WMAP 3 year data give n s = 0.951 ± 0.017 (Spergel et al. 2007). 
With these values the corrected weighted mean T is 0.209 ± 0.014. 
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Fig. 4. — The stretch parameter R and the acoustic scale l a as a function of the ratio of dark 
energy density to matter plus radiation density at zls- These curves were evaluated for flat models 
by varying w', setting wq = —0.922 — 0.309u/ and then adjusting Qm to minimize the x 2 from the 
supernovae fit plus the simple Qm = 0.296 ± 0.029 prior. 

2.6. Acoustic Oscillations 

The most important parameter from the CMB data is the acoustic scale. In this paper the 

1/2 

stretch parameter R = H Q (1 + zls)Da{zls)/ c (Bond, Efstathiou & Tegmark 1997; Wang Sz 
Mukherjee 2006) is not used, but the acoustic scale £ a evaluated at the mean baryon and dark 
matter densities is used in its place. The acoustic scale is defined as (Page et al. 2003) 

_ tt(1 + z ls )Da{zls) 
a f 1/(1+ ^c s da/(aa) ^ 

where zls is last scattering, c s is the sound speed, a is the scale factor, a is its time derivative, 
and Da is the angular size distance. The stretch parameter R has a different normalization and 

— 1/2 

approximates the denominator as oc Q, M , which is not a good approximation when the dark 
energy is significant at zls- is very well determined by the CMB data, with t a = 303.14 ± 1.04 
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Fig. 5. — Contours of \ 2 i n the ^Af-^A plane when the dark energy is assumed to be a cosmological 
constant with w = —1. The large ellipses use only the Hubble diagram data from the supernovae 
and the GRBs. The small ellipses near the "Flat" line use only the CMB acoustic scale and the 
baryon acoustic distance ratio. Contours are drawn for A% 2 = 0.1, 1, 4, and 9. The vertical lines 
show -3 . . . +3 a for the Qm prior derived from T, H and ujm- 
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Fig. 6. — Contours of x 2 m the ^m-^k plane. In the lower panel the dark energy is allowed to 
vary as w(z) = w + 2w'(l — a), using all constraints and datasets. A flat model is consistent with 
the data, and the best fit is very slight open. In the upper panel w = —1, and a flat model is again 
consistent but the best fit is slightly closed. 
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Fig. 7. — Contours of Ax 2 in the w, &k plane, with w' fixed at zero and Qm adjusted to minimize 
X 2 at each point. The black dot shows the flat ACDM model, which has 2 fewer free parameters 
and less than two units of A% 2 . 
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Fig. 8. — Contours of A% 2 in the w, w' plane, using all constraints and datasets. The dotted 
contour shows the A% 2 = 6.17 used in the DETF Figure of Merit. On the left the curvature is 
fixed at zero and £Im is adjusted to minimize x 2 at each point. On the right both Qm and are 
adjusted to minimize x 2 at each point. The black dot shows the cosmological constant w = —1. 
The dashed diagonal line shows where the dark energy is equal to the matter plus radiation density 
at last scattering. The solid diagonal line shows the 3 a limit on the stretch parameter S BBN . 

in the non-flat ACDM chain. There are correlated deviations in wg, ujm, and Q,k that produce 
correlated deviations in the numerator and denominator of Eq(|3]) which cancel out in the ratio. 
In this paper the full fit to the CMB data is not performed, but the baryon and CDM densities 
are fixed at their mean values. With this simplification, the determination of the acoustic scale 
loosens to £ a = 302.97 ±4. 14 because the deviations in the numerator are still present. The relative 
accuracy of t a with fixed wg and ujm is the same as the relative accuracy of R. The main reason 
to use £ a instead of R is that it is clear how £ a is modified when the dark energy is important at 
last scattering. Figure H] shows that t a shifts a great deal when the dark energy density is greater 
than the radiation plus matter density at recombination, while the change in R is smaller and in 
the wrong direction. 

The CMB also provides the matter density com used as part of the large scale structure analysis. 
For models far from the geometric degeneracy line, a full fit to the CMB power spectrum will 
compensate for an incorrect distance to the surface of last scattering by changing ujm and u>b to 
give a different value of the sound horizon at last scattering, but this effect has not been been 
included here. Ultimately a full fit to all the data should be performed, and the CMB is both the 
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most informative and the most computationally intensive of all the datasets. Spergel et al. (2007) 
has some limits on w, but not on the w, w' plane, and does not include either the Riess et al. (2007) 
and Wood-Vasey et al. (2007) increments to the supernova data or the GRB data. 

The second major input involving acoustic oscillations comes from the baryon acoustic oscilla- 
tions detected by Eisenstein et al. (2005). In this paper we use the ratio of the distance Dy(0.35) 
at z = 0.35 to the tangential distance at last scattering, 

^0.35 — -77—. rpr~7 r l 4 J 

(1 + zls)Da{zls) 

This ratio is easily computed even when the dark energy is significant at last scattering. The ratio 
is also slightly more precise than the A parameter since the scatter induced by the uncertainty in 
ujb and ujm cancels out in the ratio. The A parameter also uses the approximation that the sound 
travel distance is oc Qj}^ 2 which fails when dark energy dominates early. Finally, A/R is exactly 
proportional to Dy(0.35)/[(1 + zls)Da(zls)]- Using both A and the distance ratio amounts to 
double counting the baryon acoustic oscillations. 

The redshift of last scattering, zls, is essentially independent of both the baryon density and 
the expansion rate because the electron density is n e oc H(z)/a(T) while the optical depth is 
r oc <jTn e c/H(z) oc (Jt/oi{T) where <jt is the scattering cross-section and a is the recombination 
coefficient. Hence zls = 1089 is used for all calculations in this paper. 



2.7. Big Bang Nucleosynthesis 

Steigman (2006) has analyzed the light element abundances and obtained limits on a stretch 
factor S BBNS = (1 + pDE/p 1+m +k) 1/2 - The limits are S BBNS = 0.942 ± 0.030. This means that 
at 3 a the dark energy density must be less than 6.4% of the radiation density at the redshift of 
nucleosynthesis, z 10 9 . 



2.8. Expansion Rate at Last Scattering 

Zahn & Zaldarriaga (2003) have shown that even though H(z) cancels out in determining zls, 
the width in conformal time of the transition from opaque to transparent can be determined using 
the CMB TT and TE power spectra, and that this can give a stretch factor similar to the BBNS 
stretch factor. Zahn & Zaldarriaga interpret this in terms of varying Newton's constant G, but it 
really depends on \JGp. Interpreting the limit in Zahn & Zaldarriaga in terms of the dark energy 
density instead of a changing G gives S LS = (1 + Pde / P-y+m+k) 1 ^ 2 = 1-749 ± 0.471 at z LS . With 
the first year WMAP data this improved to 1.04 ± 0.24 (Zahn, private communication). 



-14- 



3. ANALYSIS 

All of the analyses in this paper depend on the expansion history and geometry of the Universe. 
The expansion history can be calculated using 



a = H y/Q. M /a + Q R / a 2 + tt K + &de[pde{z)/ Pde{0)W (5) 

The dark energy density as a function of redshift is computed using the w = Wq + w a (l — a) formula 
(Chevallier & Polarski, 2001). Following Linder (2003) I set w' = w a /2. This then gives 

Pde(z) + z)3+3wo+6w > exp / (6) 



Pde(0) \ 1 + z 

The calculation of angular size and luminosity distances then follows Wright (2006). 

Note that the binned distance moduli were not used for the analysis. Thus the Hubble constants 
used when producing the binned data tables are irrelevant in the analysis. But both the supernovae 
and the GRBs have an associated "nuisance" parameter, A4sn and Mgrb, which are adjusted to 
give the best modified x 2 at every position in the {flj,/, w, w'} parameter space. Thus 

X 2 {tt M: n k ,w,w') = min V]/((/ij - /i c (zi))/cJi) 

M * — * 

i 

with 

/i c (z) = M + 5log{D L (z;n M ,n k ,w,w')) (7) 

The form in Eq([7|) is used twice, once for the supernovae giving Xsn an( ^ once f° r tli e GRBs giving 
Xgrb • Contours of the sum of these two terms are shown in Figure [5j 

The acoustic scale quantities £ a and R0.35 give a contribution to the overall x 2 °f 

XI(^M, Mk,W, w') = f(£ a (^M, W, It/) - £ a ,obs)/cr{£a)) 

+ f(Ro.35(^M,^k,W,w') - Ro.35,obs)/cr(Ro.35), ■ (8) 

where / is the function from Eq([T]) with the transition from quadratic to linear behavior set at 3a: 
f{x) = x 2 for \x\ < 3, or 6\x\ — 9 otherwise. 

Contours of x\ are shown in Figure O 

Normally the data on T = ^m^, H q and Qm^ 2 can be combined to give a prior on £Im- There 
is an overdetermined set of equations with variables £Im and h: 




]n(h) 



( ln(0.209) ± (0.014/0.209) 

ln(0.1289) ± (0.0079/0.1289) | (9) 
^ ln(0.71) ± (0.05/0.71) 



The least squares solution of these equations gives H Q = 67.3 ± 3.7 km/sec/Mpc and FLm 
0.296 ± 0.029. Contours showing this prior are the vertical lines in Figure [5j 
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But the relationship between T = Oj^/j and the big bend in P{k) is derived assuming that 
matter and radiation are the only significant contributors to the density at z eq , the redshift of 
matter-radiation equality. Prior to z eq , the Universe is expanding faster than the free-fall time for 
matter perturbations, so growth is suppressed for fluctuations that are inside the horizon earlier 
than z eq . To allow for the possibility that dark energy contributes, z eq is found using 

n M (i + z eq f = n R (i + z eq f + n K (i + Zeq f + n DE (i + z eq f +3wo+6w ' ex P (z^i\ ( i ) 

The horizon at z eq is found using 

rl/(l+z eq ) cda 

D eq = — 11 

Jo aa 

Then the effective value of T is T e g = (1602 km / sec) / (H a D eq ) . By finding the effective F e g at 
two different values for H , which leads to two different values for Qr since £luh 2 is fixed by the 
measurement of T , the standard T = JT^f/i can be replaced by a modified power law function 
of h which reduces to the standard form except when the dark energy is significant near z eq . 
For example, in the track of models shown in Figured! when pDE/p-ym = 0.138 at zls one has 
U M = 0.2983, w = -1.098 and w' = 0.570. These values give w{z) = 0.042 at z » 10 3 , so the 
dark energy density grows faster than the matter density but slower than the radiation density at 
high z. For these parameters one finds that r e fj = 0.171704 for h = 0.675373 and T e g- = 0.188939 
for h = 0.746402. Fitting a power law function of h gives r eff m 0.2499/i 95657 instead of the 
standard 0.2983/i. For h = 0.71 this is a 16% difference which is quite significant compared to the 
7% precision of the mean of the T's from the SDSS and the 2dF. When finding x 2 , a weighted 
mean estimate for h is found using the T, H Q and co m priors. This weighted mean minimizes the 
X 2 contribution from these three priors, and this minimum is added to the x 2 from the Hubble 
diagram, CMB, and the BAO. Therefore h becomes a third nuisance parameter. Then 

Xr(£lM,tt k ,w,w') = mm[f((T cS (Q M ,^k,w,w',h) -T obs )/a(F)) 

h 

+ f(n M h 2 - u Mj0hs )/a(u:M) + f((h - h obs )/a{h))} (12) 

is added to the overall x 2 - As the dark energy becomes significant at recombination, Qm has to 
increase to keep the effective T close to the observed value. If we take the same w = —1.098 and 
w' = 0.570 example discussed above, and adjust Q,m in flat models to minimize x 2 alone, Qm 
increases by 1.5c" from 0.296 to 0.341. Thus a simple £Im prior is not correct when dark energy 
is significant at zls- Note that ujm obsi which is derived from the peak heights and trough depths 
in the CMB angular power spectrum, will be affected by perturbations in the dark energy density, 
but this effect has not been included in this paper. 

It is simple to find the BBNS S BBNS parameter for any point in the {Qm, w, w'} pa- 
rameter space, and add an appropriate term to x 2 ■ I n doing this calculation I have used h = 
iJ(jJm{CMB)/£Im for simplicity instead of combining this analysis with the calculation of x 2 - 
Since the desired value for S BBN is slightly less than one, while dark energy can only increase 
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S BBNS above one, the BBNS data act only as an upper limit on w' . The Zahn & Zaldarriaga limit 
on dark energy at recombination is handled the same way. This gives a stretch factor term of 

,,2 _ ffQBBNS(C) O .,A qBBNS\ ii qBBNS\\ 
Xs - J\ b {^hl,^k,w,w ) - b obs )/cr{b >) 

+ f(S LS (n M ,n k ,w,w') - S^ b s s )/a(S LS )). (13) 

The final form for the overall x 2 is 

X 2 (0 M , fijfe, w, w') = Xsn + Xgrb + xl + Xr + x|- (14) 

Note that a minimization over Msn has already been done in computing Xsn ( see ^l 11 EJ) j a 
minimization over Mgrb has already been done in computing Xgrb > an d a minimization over /i 
has already been done in computing xp (see Eqn ll2p . Thus the overall x 2 is 

X 2 (tt M ,tt k ,w,w') = min x 2 (^M,tt k ,w,w' \h, M S n , M G rb) (15) 

h,MsN,MGRB 

All of the contour plots in this paper except for Figure [5] show contours of this combined function. 
Unplotted variables are either fixed at assumed values or minimized over. It is correct to remove 
nuisance parameters by minimizing x 2 or maximizing the likelihood over them (Cash 1976). 

Marginalization by integrating the likelihood over the nuisance parameters is wrong (Wright 
1994), although it is allowable to marginalize by integrating over the a posteriori probability density 
function. But this requires a correct prior distribution for the nuisance parameters. The Monte 
Carlo Markov chain (MCMC) technique is a useful tool for performing such integrals, but the 
MCMC thus requires correct priors. As an example, consider Figure 20 of Spergel et al. (2007) 
which shows the CMB constraint in the (Omj^a) plane. This should have a uniform prior in 
(fljv/, ^a)> but the chain was computed with a uniform prior in Q s ), so the chain weights had 
to be multiplied by the Jacobian of the transformation between (f2A,0 s ) and ((]m,^a) before a 
subset of the chain was selected for plotting. 

Marginalization by minimization avoids the need for a prior on the unplotted variables. The 
plots in this paper show contours of the likelihood, which are the same as contours of the a posteriori 
probability density function if one assumes a uniform prior in the plotted variables. 

4. DISCUSSION 

There are a large number of different cuts through the 4 dimensional parameter space that can 
be plotted, and when different subsets of the data are considered the number of plots multiplies 
rapidly. The best fit x 2 f° r some of these combinations are listed in Table HI Figure [5] shows the £Im 
vs. Q\ plane when the dark energy is constrained to be a cosmological constant with w = — 1. Three 
data subsets are shown: the supernova and GRB Hubble diagrams, the CMB plus BAO acoustic 
scale data, and the T, H Q and u>m data. Clearly the acoustic scale data give a strong confirmation 
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Table 4: x 2 for fitting to 358 SNe, 69 GRBs, T, H a , uj m , and BBNS. The parameters in the fit are 
w, w', U M , 9, K , M S n, M G rb, and h. 



Model type 


w 


w' 




n K 


x 2 


flat ACDM 


-1 





0.306 





427.905 


nonflat ACDM 


-1 





0.315 


-0.011 


426.983 


nonflat constant w 


-0.894 





0.309 


-0.003 


426.249 


flat varying w 


-1.126 


0.451 


0.305 





423.580 


nonflat varying w 


-1.098 


0.506 


0.299 


+0.015 


422.856 



of the need for dark energy, and a much better constraint on the curvature of the Universe. Using 
all the data together gives the plot shown in Figure El The best fit model is slightly closed with 
&tot = 1-011 and VLm = 0.315. The best fit flat ACDM model has less than one more unit of x 2 
than the best fit non-flat ACDM model so there is no evidence for spatial curvature from these fits. 
Figure [U] also shows the effect of allowing w(z) to vary. 

Another way to see this is to plot \ 2 vs - w an< l ^if> as seen in Figure [71 In this plot w' is 
forced to zero, and Qm is adjusted to minimize \ 2 a t each point. This plot is very similar to a 
comparable plot in Spergel et al. (2007). 

The flat ACDM model has only 5 more units of \ 2 than a non-flat variable w model with 
3 more free parameters. The probability of this occurring by chance alone is over 16%, so this 
improvement is not significant. 

Plots of the w vs. w' plane are shown in Figure EJ with and without the assumption of a flat 
Universe. It is obvious that the contours would have extended to much higher values of w' if the 
nucleosynthesis and acoustic scale constraints had not been used. The tilts of the ellipses below 
these cutoffs indicate that the pivot redshifts for the priors and datasets used here are z = 0.4 and 
z = 0.22. Values of w at the pivot redshift are close to —0.9 which is the best fit when w' is forced 
to zero, as shown in Figure [71 

The Dark Energy Task Force (Albrecht et al. 2006) has defined a Figure of Merit as one over 
the area of the A% 2 = 6.17 contour in the w,w a plane. Albrecht et al. said that this contour is 
a 95% confidence contour, which not exactly correct. It is actually a 95.4% confidence contour to 
match the 2<r confidence of a one dimensional Gaussian (Albrecht, private communication). The 
FoM for the flat case is 1.59 while for the non-flat case it is 1.39, but both of these cases have run 
up against the constraints against dark energy domination before last scattering. These constraints 
have given the non-flat case a bigger boost than the flat case. In fact the area of the A% 2 = 1 
contour is smaller for the non-flat case than the flat case, which shows both the inability of the 
current data to constrain w, w a and the importance of using the constraints coming from high 
redshift physics. 
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5. CONCLUSION 



All of the data considered here are consistent with a flat ACDM model with a constant equation 
of state w = —1. The current data, even with well over 300 supernovae, are not adequate for 
measuring a time variable equation of state with reasonable precision. Serra, Heavens & Melchiorri 
(2007) and Davis et al. (2007) agree with this conclusion, as do Liddle et al. (2006), Alam et al. 
(2006) and Li et al. (2006) using earlier and smaller datasets. The current acoustic scale data, seen 
in the CMB and the baryon oscillations, is giving more precise information about the expansion 
history of the Universe, but without the dense redshift coverage provided by the supernovae. There 
appears to be a systematic deviation of the supernovae data from the models around redshifts near 
0.5, whose origin is unknown. Since the choice of data subsets affects the size of this deviation it is 
probably an artifact. Nesseris & Perivolaropoulos (2006) have found systematic differences between 
the different data sources that went into Riess et al. (2007) dataset, so artifacts are not unlikely. 
Furthermore the scatter between the distance determinations of identical supernovae by different 
groups is unexpectedly large. The GRB Hubble diagram is not very precise but it appears to be 
consistent with the supernova data. The GRB Hubble diagram does help break the degeneracy 
between w ^ — 1 and Qk ^ 0. The 2dF and SDSS values for T differ by a slightly disturbing 
amount, so a new and improved measurement of T would be useful as a tie breaker. The Hubble 
constant appears to be determined to better than 10% but independent new data would be quite 
valuable in pinning down 0,m in combination with T. Better CMB data from Planck should reduce 
the uncertainty in and wm, which will reduce the uncertainty in R or l a . It is clear that better 
data of many types will be needed to pin down w(z). 
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